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Brooks and Makover introduced an approach to studying the 
global geometric quantities (in particular, the first eigenvalue of the 
Laplacian, injectivity radius and diameter) of a "typical" compact 
Riemann surface of large genus based on compactifying finite-area 
Riemann surfaces associated with random cubic graphs; by a theo- 
rem of Belyi, these are "dense" in the space of compact Riemann 
surfaces. The question as to how these surfaces are distributed in the 
Teichmuller spaces depends on the study of oriented cycles in random 
cubic graphs with random orientation; Brooks and Makover conjec- 
tured that asymptotically normalized cycle lengths follow Poisson- 
Dirichlet distribution. We present a proof of this conjecture using 
representation theory of the symmetric group. 

1. Introduction. The study of the first eigenvalue of the Laplace operator 
on compact Riemann surfaces of increasing genus has received considerable 
attention over the last thirty years; see [18] and references therein. On the 
one hand, we have a celebrated theorem of Selberg [48] (see [36] and [32] 
for refined estimates toward Selberg's conjecture) and its generalization by 
Sarnak and Xue [47] (see [29] and [26] for related results), asserting that 
the first eigenvalue of the congruence surfaces of arbitrary genus is bounded 
away from zero; on the other hand there are examples due to Selberg [48], 
Randol [44] and Buser [16], showing that, in general, the first eigenvalue can 
be made arbitrarily small. The surfaces in the latter examples are "long and 
thin," so in some sense live on the boundary of Teichmuller spaces, and it is 
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a fascinating question to determine what happens for a "typical" Riemann 
surface of large genus. 

The idea of using cubic graphs to study the first eigenvalue of Riemann 
surfaces originated in the work of Buser [15, 17]. As we discuss in Section 3, 
the behavior of the first eigenvalue of the discrete Laplacian on a random cu- 
bic graph is understood rather well. In [14] , Brooks and Makover introduced 
an approach to studying the first eigenvalue of the Laplacian of a "typical" 
compact Riemann surface of large genus based on compactifying finite-area 
Riemann surfaces associated with random cubic graphs with random orien- 
tation; we review the Brooks-Makover construction in Section 2. To each 
cubic graph T with an orientation O, they associate two Riemann surfaces: 
S (r, O), a finite-area noncompact surface, and S c (T,0), a compact sur- 
face. The surface S°(T,0) is an orbifold cover of H/PSX(2,Z) described 
by (r,C) and therefore shares some of the global geometric properties with 
the graph T. The compact surface S (T,0) is a conformal compactification 
of S°(T,0); Brooks and Makover proved that almost always the global ge- 
ometry of S (r,C) is controlled by the geometry of S°(T,0). Moreover, 
according to Belyi's theorem [4], the surfaces S C (T, O) are precisely the Rie- 
mann surfaces which can be defined over some number field and so form a 
"dense" set in the space of all Riemann surfaces. 

The question as to how these surfaces are distributed in the Teichmiiller 
spaces depends on the study of oriented cycles in random cubic graphs 
with random orientation; Brooks and Makover conjectured that asymptoti- 
cally normalized cycle lengths follow Poisson-Dirichlet distribution. We re- 
call the definition of Poisson-Dirichlet distribution [3]. Let Bi,B 2 ,... be 
independent random variables uniformly distributed on [0, 1] . Define G = 
(Gi, G 2 , • • •) as follows: 

G 1 =B 1 ; G 2 = (1-B 1 )B 2 ; G 4 = (1 - B x )(l - B 2 ) ■ ■ • (1 - Si-i)^. 

The random sequence G can be viewed as a description of a random breaking 
of a stick of unit length into an infinite sequence of subintervals. A stick of 
length B\ is broken off at the left, which leaves a piece of length 1 — B\. 
From this, a piece of length (1 — B\)(l — B 2 ) is broken off, and so on. G is 
a distribution on the set 



n 



x 6 R°° : x\,x 2 , • • • > 0, Xi = l 



=1 



The ranked version of G, (G^ , G( 2 ) )•••)> wnere Gm > G7 2 ) > • ■ • , has Poisson- 
Dirichlet distribution. If, instead of uniform distribution, Bi has beta-(l,0) 
density 9(1 — x) e ~ l on [0, 1] with 9 > 0, the resulting distribution is called 
Poisson-Dirichlet distribution with parameter 9. Poisson-Dirichlet distribu- 
tion arises in a great variety of problems; see [3, 42] and references therein. In 
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a recent breakthrough work [21], Diaconis, Mayer- Wolf, Zeitouni and Zerner 
proved a conjecture of Vershik [52] asserting that Poisson-Dirichlet distri- 
bution is the unique invariant distribution for uniform split-merge transfor- 
mations. 

As we discuss in Section 3, the distribution of cycle lengths for random 
regular graphs in any of the standard models is a rather difficult problem; 
adding random orientation, at first, seems only to complicate matters (see 
the example in Section 2) . However, as we explain at the end of Section 3, this 
extra randomness turns out to help us: by giving a permutational model for a 
random regular graph with random orientation, we can convert the problem 
into one involving the symmetric group. More precisely, the distribution 
of oriented cycles in a random /c-regular graph on n vertices with random 
orientation is the same as the distribution of cycles in the permutation fia, 
where [3 is chosen uniformly on the conjugacy class consisting of the product 
of A:-cycles, and a is chosen with uniform probability on the conjugacy class 
consisting of the product of 2-cycles in the symmetric group Sn with N = 
nk. In Section 4, using the Diaconis-Shahshahani upper bound lemma [22], 
the estimate on the number of rim hook tableaux by Fomin and Lulov [23] 
and representation theory of the symmetric group (in particular, the hook 
length formula and the Murnaghan-Nakayama rule) , we show that as n — > oo 
the distribution of (3a converges to uniform distribution. We then invoke 
what is perhaps the oldest occurrence of Poisson-Dirichlet distribution — 
the distribution of normalized cycle lengths for a random permutation in S n 
as n tends to infinity [49, 57, 58] — to prove the conjecture of Brooks and 
Makover. 

It turns out that the number of oriented cycles in random cubic graphs 
with random orientation was also studied by Pippenger and Schleich [43] 
in connection with topological characteristics of random surfaces generated 
by cubic interactions. The surfaces considered by Pippenger and Schleich 
are obtained by taking 3n arcs of an even number of oriented triangles 
and randomly identifying them in pairs respecting the orientation; these 
surfaces arise in various contexts in two-dimensional quantum gravity and as 
world sheets in string theory. Random cubic graphs with random orientation 
provide an alternative way of constructing the surfaces in [43] . Denoting by 
/ (n) the number of oriented cycles in a random cubic graph on n vertices with 
a random orientation, Pippenger and Schleich proved that E(Z(n)) = logn + 
O(l) and Var(Z(n)) = O(logn). Further, based on empirical study of 10,000 
random surfaces, each constructed from 80,000 triangles, they conjectured 
that 

(1.1) E(Z(n))=log(3n) + 7 + 0(1) 

and 



(1.2) 



Var (Z(n)) = log(3n) + 7 - vr 2 /6 + o(l), 
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where 7 = 0.5772 ... is Euler's constant. It was pointed out by Nicholas 
Pippenger that (1.1) and (1.2) follow from the main theorem of this paper 
and a priori bounds on the first few moments of l(n) obtained by the methods 
of [43]. Another consequence of the main theorem (Corollary 5.2) is that 
the expected area of the largest embedded ball in a random Belyi surface 
converges to of the total surface area. We conclude Section 5 by briefly 
discussing the conjectured Tracy- Widom distribution for the second largest 
eigenvalue of random regular graphs; we hope that our results will be useful 
in approaching this fascinating open problem. 

2. Belyi surfaces. In [4], Belyi proved a remarkable result asserting that a 
Riemann surface S can be defined over the field of algebraic numbers Q if and 
only if there exists a covering / : S — > C unramified outside {0, 1, 00}. We call 
such surfaces Belyi surfaces. In this section we review the Brooks-Makover 
construction of Belyi surfaces from cubic graphs. We remark that Mulase 
and Penakava [38] have given a very interesting alternative construction of 
Belyi surfaces; in their construction, the edges of the graphs are allowed to 
have variable lengths. 

Let r be a cubic graph. An orientation O on the graph V is an assignment 
for each vertex v of T of a cyclic ordering of the half-edges incident to 
that vertex. Given a cubic graph on n vertices, it will have 2 n different 
orientations. A left-hand-turn path (LHT path) on V is a closed path on T 
such that, at each vertex, the path turns left in the orientation O. 

Given a pair (T,0), we construct a finite- area Riemann surface S°(T, O) 
as follows. We take the ideal hyperbolic triangle T with vertices 0, 1 and 




Fig. 1. The marked ideal triangle T. 
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oo shown in Figure 1. The solid lines in Figure 1 are geodesies joining the 
points i,i + 1 and with the point 1+t ^ , while the dotted lines are horo- 
cycles joining pairs of points from the set {i,i + 1, ^p}. We may think of 
these points as "midpoints" of the corresponding sides of the ideal triangles, 
even though the sides are of infinite length. We may also think of the three 
solid lines as edges of a graph emanating from a vertex. We may then give 
them the cyclic ordering + ^jp). Given a cubic graph T with an ori- 
entation O, we construct S°(T,0) by associating to each vertex an ideal 
triangle, and gluing neighboring triangles. We glue two copies of T along 
the corresponding sides, subject to the following two conditions: 

(a) the "midpoints" of the two sides are glued together, 
and 

(b) the gluing preserves the orientation of the two copies of T. 

The conditions (a) and (b) determine the gluing uniquely. It is easily seen 
that the surface S°(T,0) is a complete Riemann surface with a finite area 
equal to irn, where n is the number of vertices of T. 

It is easy to see that the horocycle pieces on each T glue together to 
give closed horocycles about a cusp of S (T, O). The length of each closed 
horocycle is precisely the length of corresponding LHT path. The surface 
S (r,C) is the conformal compactification of S°(T,0). The oriented graph 
(r,C) describes S°(T,0) as a covering space of HI 2 /.P,5X2(Z), with each 
vertex being an orbifold point of order 3. 

We are now ready to state the following result, whose proof was sketched 
by Brooks and Makover in [14]: 

Lemma 2.1. S is a Belyi surface if and only if S = S (V, O) for some 
cubic graph V . 

Proof. We first show that if G is a torsion- free finite index subgroup 
of PSL2CZ1), then M 2 /G = S (T,0) for some cubic graph T. Indeed, we can 
take as a fundamental domain for PSL<i(7fc) a set F = {0 < < 1, \z\ > 
1, \z— 1| > 1}. Three copies of F fit together to give a marked ideal triangle as 
presented in Figure 1; they are transformed by means of an elliptic element 
w of order 3 in PSL2 (Z) . 

Now since G is torsion free, F, w(F) and w 2 (F) are not equivalent under 
G, and so they can all be included in the fundamental domain of G. In 
particular, there is a fundamental domain of G such that consisting of the 
copies of ideal triangle T. The graph dual to this triangulation, together with 
the boundary pairings and orientation of M 2 /G, is exactly the pair (T,0). 

Now since S is a Belyi surface if and only if one can find finitely many 
points {pi, . . . ,pi} on S such that S — {p±, . . . ,pi} is isomorphic to M 2 /G, 
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where G is a torsion- free finite index subgroup of PSL^i^) [31], the lemma 
is proved. □ 

We define probability on the space of oriented graphs with n-vertices 
(T n ,0) as follows. We pick a random cubic graph with n vertices using the 
Bollobas model, described in the next section, then we pick an orientation O 
with equal probability from all 2 n possible orientations on the given graph. 
If Q is a property of graphs, we denote by Prob n [(5] the probability that 
an oriented graph (T n ,0) picked from our probability space has property 
Q. Brooks and Makover proved [14] the following result on Belyi surfaces 
constructed from random cubic graphs with random orientation: 

Theorem 2.1. There exist positive constants C\, Ci, C3 and C4 such 
that, as n — > 00: 

(a) The first eigenvalue \\(S C (T,0)) satisfies 

ProMAi(S c (r,0))>(7i]->l. 

(b) The Cheeger constant h(S c (T,0)) satisfies 

Prob n [h(S c (T,0))>C 2 ]^l. 

(c) The shortest geodesic syst(S c (T, O)) satisfies 

Pr 6 n [syst(S c (r,O))>C 3 ]^l. 

(d) The diameter diam(S ' (r , O)) satisfies 

Prob n [diam(S c {T,0)) < C 4 log (genus (5 c (r, O)))} -» 1. 

The theorem is proved, roughly, as follows. First, one establishes that, 
with probability tending to 1 as n — > 00, a property in question holds for 
random cubic graphs [see Section 3 for results pertaining to parts (a) and 
(b)]. Then, using the fact that (T,0) describes S°(F,0) as an orbifold cov- 
ering, one transfers this information to open surfaces S°(T,0), using the 
results of Brooks in [8, 9]. One then transfers the desired property to the sur- 
faces S c (T,0) by using the Ahlfors-Schwarz-Wolpert lemma as developed 
by Brooks in [10, 11] and extended by Brooks and Makover in [12, 13, 14]. 

The topology of the surface can be read off from (T, O), using LHT paths. 
In particular, the genus is given by Euler's formula; as is customary when 
using Euler's formula, we refer to an LHT path as a face of the oriented 
graph. 

Let n be the number of vertices of T and let l(T,0) be the number of 
disjoint left-hand-turn paths. For a cubic graph, the number of edges is 



POISSON-DIRICHLET DISTRIBUTION FOR RANDOM BELYI SURFACES 7 




Fig. 2. Changing orientation on the cube. 



£(T) = ^ and the number of faces is T = l(T,0). Therefore, by Euler's 
formula, we have: 

n — 21 

(2.1) genus(5 c (r n , O)) = genus(5°(r n , O)) = 1 + — ^— . 

According to Euler's formula, in order to find the genus of a surface con- 
structed from a cubic graph with n-vertices, we need to estimate the number 
of faces (i.e., left-hand-turn paths). As explained in [14], the length of the 
largest LHT path L also determines the area of the largest embedded ball. 

As we discuss in Section 3, computing the distribution of closed paths in 
random regular graphs is a very difficult problem. An additional complica- 
tion is that left-hand-turn paths are not necessarily simple closed paths on 
the graph, but can self-intersect in a very complex pattern. We would like to 
illustrate these complications with the simple example of the 1-skeleton of 
the cube, with the usual orientation (Figure 2A). We have J- =6, and all the 
faces are simple paths of the graph, therefore the genus is 0. In Figure 2B 
we changed the orientation of the upper right vertex. With the new orienta- 
tion, the three simple faces that were previously adjacent to the upper right 
vertex are now joined to one composite face, while the other three faces are 
unchanged, hence J- =4, and the surface has genus 1. 

3. Random regular graphs. In this section we briefly review the per- 
tinent facts on random /c-regular graphs; see Wormald's paper [59] for an 
excellent survey. Given a /c-regular graph Q and a subset X of V, the expan- 
sion of X, c(X), is defined to be the ratio |<9(X)|/|X|, where d(X) = {y £ 
Q : distance^, X) = 1}. The expansion coefficient of a graph Q is an analogue 
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of Cheeger's constant for Riemann surfaces and is denned as follows: 

(3.1) c(g) = M{ c (x)\\x\<±\g\}. 

A family of /c-regular graphs X n ^ forms a family of C-expanders [34, 46] if 
there is a fixed positive constant C such that 

(3.2) liminfc(X nfe ) > C. 

The adjacency matrix of Q, A(Q), is the \G\ by \G\ matrix, with rows and 
columns indexed by vertices of Q, such that the x,y entry is 1 if and only if 
x and y are adjacent, and is otherwise. For a fc-regular graph, the adjacency 
matrix is related to the combinatorial Laplacian A by A = kl — A. Using the 
discrete Cheeger-Buser inequality, condition (5.1) can be rewritten in terms 
of the second largest eigenvalue of the adjacency matrix A{Q) as follows: 

(3.3) limsup A2(Ai,fc) < k. 

n— >oo 

In 1973, Pinsker [41] observed that a random regular graph is a good 
expander. This corresponds to the following fact about random matrices: a 
random symmetric matrix of size N with k ones in each row and column 
and all other entries zero has its biggest eigenvalue equal to k, but its next 
eigenvalue will be bounded away from A; by a fixed amount independent of N. 
The result of Pinsker on expansion coefficients of random regular graphs was 
considerably strengthened by Bollobas [7] , who also introduced a widely used 
configuration model for random regular graphs [6]. In this model, random 
/c-regular graphs on TV" vertices are represented as the images of so-called 
configurations. Let W = Uj=i Wj be a fixed set of 2m = nd vertices, where 
\Wj\ = d. A configuration F is a partition of W into m pairs of vertices 
called edges of F. Clearly there are 

(3.4) N=N(m)= (^)(^)...(l)/ m , = ^ 



2 J V 2 J \2 r 

configurations. [We write (a)b = a(a — 1) • • • (a — b + 1).] 

Let $ be a set of configurations; we turn it into a probability space by 
giving all configurations the same probability. We now define a map $ — > Q n k 
as follows. Given a configuration F, let 4>{F) be the graph with vertex set 
V = 1, . . . ,n in which ij is an edge iff F has a pair with one end in Wi 
and the other in Wj. Every G £ G n ^d is the image of G = 4>(F) for (dl) n 
configurations. The number of configurations containing a given fixed set of 
I edges is 

iV,W=( 2 " , - 2, )( 2 " 1 - 2 2, - 2 )-(^)/( m -/)! 

N(m) 



(2m - l)(2m - 3) ■ • ■ (2m -21 + 1) 
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Using the configuration model, and in particular (3.5), Bollobas proved 
the following result: 

Theorem 3.1 ([5]). Let Xi denote the number of closed walks of length 
i in Q G Gnk- Then, for i = o(logfc_ 1 ?i) as n — > oo, random variables Xi 

(k—lY 

converge to independent Poisson random variables with mean v 2 . . 

Counting cycles of length greater than logn is substantially more diffi- 
cult. In a recent breakthrough work [25], following his earlier work in [24], 
Friedman estimates the number of cycles of length 0(log 2 n) and uses this 
estimate (among other things) to prove that /c-regular graphs on n vertices 
G n ,k are asymptotically Ramanujan: for k fixed and e > 0, the probability 
that Ai(X nj fc) < 2\/k — 1 + e tends to 1 as n — > oo. 

The bound of 2y/k — 1 is optimal in view of the result of Alon-Boppana 
[1, 35]. We also mention an early result of McKay [37], who showed that 
spectral density of random /c-regular graphs converges to Kesten's measure, 
that is, a measure supported on [— 2\/k — 1, 2y/k — 1 ] and given by 

k ^/4(k-l)-t 2 



(3.6) v k 



2tt k 2 - t 2 



We now introduce the permutational model for random regular graphs 
with random orientation. Consider a neighborhood of a vertex of a cubic 
graph. It contains three half-edges incident to it. Denote the set of all half- 
edges by H; \H\ = 2\E\ = 3|V|. Now the cyclic ordering of half-edges at each 
vertex is specified by a 3-cycle, and a collection of all cyclic orders on the 
half-edges yields a permutation of H consisting of the product of 3-cycles. 
The structure of the underlying graph is given by the way half-edges couple 
to each other to form an edge; this is described by a permutation all of whose 
cycles are of order 2. 

In the Figure 3 we have 



/3 = (1,3,5)(2,12,8)(4,7,9)(6,10,11), 
a =(1,2) (3, 4) (5, 6) (7, 8) (9, 10) (11, 12) 

and 

(p = /3 a = (1, 12, 6)(2, 3, 7)(4, 5, 10)(8, 9, 11). 

Let Q G G n ,k be a fc-regular graph on n vertices; let O be its orientation, 
that is, a cyclic ordering of incident half-edges for each vertex. Associated 
with a (G, O), we have a pair of permutations (f3,a) in Sn, where N = kn, 
with cycles of (3 encoding the information about the vertices with their cyclic 
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orientation and a encoding the information about the edges. The cycles of 
<p = (3a encode the information about the faces. 

Now the problem of the distribution of the lengths of LHT paths can be 
given the following, equivalent, formulation. Denote by C r the conjugacy 
class of A jv consisting of the product of N/r disjoint r-cycles. Choose (3 
with uniform probability on Ck and choose a with uniform probability on 
C2. Then the distribution of LHT paths in a random /c-regular graph on n 
vertices with random orientation is the same as the distribution of cycles in 
the permutation (3a. 

4. Proof of Brooks Makover conjecture. This section is devoted to the 
proof of the following result: 

Theorem 4.1. Let N = kn and let P r denote the probability measure on 
An supported on C r ; let U denote the uniform distribution on An- Then, 
for k>3, 

(4.1) \\P k *P 2 -U\\ — 0. 

n^oo 

Here 

(4-2) \\f-g\\=wsx\f(A)-g{A)\ 
is a total variation distance. 
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Corollary 4.1. The distribution of LHT paths for random regular 
graphs with random orientation converges to Poisson-Dirichlet distribution; 
in particular, the Brooks-Makover conjecture is true. 

Proof. Let it denote the cycle type of a permutation; denote by C s (ir) 
the number of permutations with a given cycle type tt in the symmetric 
group Sn and by C a {tt) the number of permutations with cycle type tt in 
the alternating group A^. As is well known, 

C s (7r) = C s ( ai ,...,a N ) 

(4.3) 



(4.4) C a (tt) 



I], i^Wa^vr)! 

for 7r = l ai ... N aN . The cycle index polynomial of Sat is given by 
p N ( Xl ,...,x N )= ]T C s (ai,...,a N )x ai ---x aN . 

laiH \-Na N =N 

It is easy to see that the cycle index polynomial of An, qn{^u ■ ■ ■ > x n), is 
given by 

q N (xi, . . .,x N ) =p N (x 1 ,.. .,x N ) +pn(x 1 ,-x 2 , • • • , (-l)^ -1 ^). 
It follows that 

2C s (tt), if a2 + • • • + aw is even, 
0, otherwise. 

Now the condition that 02 + • — h is even, is exactly the condition for 
a permutation to be in An and, subject to this condition, we see that the 
distribution of cycles is given by exactly the same formula as in the case of 
the symmetric group. Consequently, the proofs of the asymptotic Poisson- 
Dirichlet distribution for normalized cycle lengths in the symmetric group 
apply also in the case of alternating group. The proof of the corollary is 
completed by applying Theorem 4.1 and the triangle inequality. □ 

We now turn to the proof of Theorem 4.1. The basic tool is the following 
result, known as the Diaconis-Shahshahani upper bound lemma; see [19, 20] 
for a survey of applications and ramifications: 

Proposition 4.1 ([22]). Let G be a finite group and denote by G the set 
of irreducible unitary representations of G. Let P be a probability measure 
on G and denote by P{p) its Fourier transform at the representation p£ G. 
Let U be uniform probability measure on G. Then 



(4.5) ||P - U\\ 2 < \ ]T dim(p)tr(P( /0 )P(p) 



peG 
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SYT T 



Here ||P — U\\ is the total variation distance denned in (4.2) and 'id' 
denotes the trivial representation. 

We apply Proposition 4.1 to the situation where G = An and P = P k * P 2 . 
Since Pk is a class function, Pk(p) is a scalar matrix given by X p (CfcKdim(p)) 
where x p (Cfc) is a character of the alternating group associated with repre- 
sentation p evaluated on the conjugacy class and idim(p) is the identity 
matrix of dimension dim(p). Now the Fourier transform maps the convolu- 
tion of functions Pk*P 2 into their product Pk * P2 = PkP2 and, consequently, 
(4.5) implies 



(4.6) 



\P k *Pi-U\ 



1 

< - 
- 4 



E 



X p (C k ) X p (C 2 ] 
dim(p) 



The representation theory of the alternating group An is closely allied 
with the representation theory of the symmetric group Sn [30] . Represen- 
tations of the symmetric group Sn are labeled by partitions A h N. A parti- 
tion A of a nonnegative integer N is a sequence (Ax, • • • , A r ) G N r satisfying 
Ax > • • • > A r and A* = iV. We call |A| =J2^i * ne s ^ ze °f ^- The number 
of parts of A is the length of A, denoted /(A). We write rrii = rrii(X) for the 
number of parts of A that are equal to i, so we have A = (l m i2 m2 . . .). 

The Young diagram of a partition A is defined as the set of points (i,j) G 
1? such that 1 < i < Xj] it is often convenient to replace the set of points 
above by squares. The conjugate partition A' of A is defined by the condition 
that the Young diagram of A' is the transpose of the Young diagram of A; 
equivalently mj(A') = A, — Aj+x- 

A standard Young tableau (SYT) of shape A is a filling of the boxes of 
A with positive integers 1, ...,N such that the rows and the columns are 
strictly increasing. 

In Figure 4 we exhibit a partition A = (5,5,3,2) = (1°2 1 3 1 5 2 ) and an 
SYT T of shape A. 

We denote by f x the dimension of representation of Sn indexed by A; it is 
equal to the number of standard Young tableaux of shape A [30]. We denote 
by x A (aO t ne value of the character indexed by A on the conjugacy class 
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Cfj,. For example, for the conjugacy class Cfe, the corresponding partition is 
H = k n . 

The partition A is called self-associated if its Young diagram is symmetric 
with respect to the main diagonal, equivalently if A = A'. Now for an irre- 
ducible representation of SV indexed by A which is not self-associated, the 
restriction of A to Ajy is irreducible ([30], page 67). If A is self-associated 
then restriction to An splits into two irreducible representations, A + and 
A~ . The character values of the irreducible representations A + and A - can 
be expressed in terms of character values of x A (m); f° r ^ ne conjugacy class \i 
not equal to the set of main diagonal hooks of A (see the definition of hook 

below), it is given by x ^ . All irreducible representations of are ob- 
tained from the irreducible representations of Sn in this way. Consequently, 
we obtain the following lemma: 



Lemma 4.1. Notation being as above, we have 

(«) £(™^ E (*°>W' 



dim(p) / \ y 

p&A N AeSjv 



Proof. Indeed, if A is not a self-associated partition, we have dim J 4 iV (A) = 
/ A and x\ N ip) = X A (/-*)• If A is a self-associated partition, then, since con- 
jugacy class Ck is not an array of main diagonal hooks of A, we have 

dim J 4 JV (A ± ) = 4p and XAwi^k) = x ^ ■ This completes the proof of Lemma 
4.1. □ 



We thus have to analyze the sum 

AhTV \ J / 

We briefly recall the basic facts pertaining to f x and x A (m)j referring to 
[30] for more details. Given a diagram A and a square u = G A, a hook 
with vertex u is a set of squares in A directly to the right or directly below 
u. We define hook length (also referred to as hook number) h{u) of A at u by 

h{u) = Xi + X'j — i — j — 1. 

Equivalently, h(u) is the number of squares directly to the right or directly 
below u, counting u itself once. For instance, in Figure 5 we display hook 
lengths for the partition A = (5, 5, 3, 2). 



14 



A. GAMBURD 



Theorem 4.2 (Hook length formula). Notation being as above, we have 

, N\ 
(4.9) / - 



Projecting a hook with vertex u onto the boundary (rim) of A yields a 
rim hook. By definition, a rim hook is a connected skew shape with no 2 x 2 
square. The height of a rim hook, ht(-R), is defined to be one less than its 
number of rows. 

A rim hook tableau T of shape A and type //, where fx = (/ii, . . . ,fJ, n ), is 
an assignment of positive integers to the squares of A such that every row 
and column is weakly increasing; the integer i appears /ij times, and the set 
of squares occupied by i forms a rim hook. The height of a rim hook tableau 
is defined to be the sum of the heights of rim hooks appearing in T. 

In Figure 6, we exhibit a rim hook tableaux of shape A = (5,5,3,2) and 
type fi= (4,4,4,3). 

In the particular case when all the parts of \x have size k, the tableaux 
described above is referred to as a k-rim hook tableaux] we denote the num- 
ber of /c-rim hook tableaux by ft. In particular, f± = f x is the number of a 
standard Young tableaux. 

Theorem 4.3 (Murnaghan-Nakayma rule). Notation being as above, 
we have 



^) = E(-i)^ (T) > 

T 

where the sum is over all rim hook tableaux of shape A and type \x. 



(440) **(/*) =£(-1)' 

T 
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Fig. 5. 
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For fx = k n , all signs (4.10) are the same [30], page 82; consequently, 

(4-11) \ X \C k )\ = fl 

In the proof of Theorem 4.1 we will make crucial use of the following 
estimate for f£: 

Theorem 4.4 ([23]). Suppose N = kn. Then 

(4.12) ft < (£^(/ A ) lA = 0(^-V(^) )(/ A ) iA ) 

where the implied constant depends only on k and not on A. 

Another crucial ingredient in the proof of Theorem 4.1 is furnished by the 
following proposition: 

Proposition 4.2. For any t>0, we have 

(4.13) £ (/*)"* = 0(tf _mt ), 

XhN 

where the implied constant depends only on m. 

We remark that the sums of the form J^w-N^f^) 13 f° r > have been 
studied by Regev [45] and Vershik [56]. In particular, Regev relates the 
asymptotic computations of such sums to the matrix integral of the form 

/oo poo n 

••• / J\x i -x J fe-W 2 ^+- +x N)dx 1 ---dx N , 
-oo J-oo i<: j 

this being one of the first hints of the deep connection between random 
matrices and random permutations; see [40] and references therein for a 
recent survey. 

Proof of Proposition 4.2. First, we observe that since f x = f x ', it 
suffices to prove Proposition 4.2 for the sum 

Ah TV 
A'^Ai^JV-m 



Now we split this sum into three parts: 
(4-14) Ei= E (/ 



A\-t 



AhiV 

a;<Ai 

3A r /4<Ai<AT-m 
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(4-15) S 2 = £ (/V 

XhN 

7V/8<Ai<37V/4 

and 

(4-16) S 3 = Yl 

AhTV 
A'j<Ai 

Ai<jV/8 

To analyze Si, we first note that 
(4.1T) /^(A) ifAl> f' 

Indeed, since the product of hook numbers for a partition of size Z is at 
most U, the product of hook numbers in the rows below the top one is at 
most (N — Ai)!. The set of hook numbers in the top row consists of distinct 
numbers not exceeding N. Since the length of the second row is at most 
r = N — Ai, the hook numbers in the rightmost \\ — r cells of the first row 
are 1, . . . , 2Ai — N . The product of the remaining r hook numbers in the first 
row is no greater than N(N — 1) . . . (N — r + 1) = j^. Applying the hook 
length formula (4.9) then completes the proof of the estimate (4.17). 
Continuing to denote r = N — Ai, we can now estimate Si using (4.17): 

m<r<N/A { r ) 

where p(r) is the number of partitions of r. Since the number of partitions 
p(r) satisfies the following inequality [2] valid for all r > 1: 



(4.19) p(r) < exp(7iy 2r/3 ), 

we have 

m<r<7V/4 I r ) 



for absolute constant c\ = e^v 2 / 3 . Let a r = ^ N \ r ^ ■ We have 

(N-r)(r + l) 
a ' r+1 ~ {N-2r){N-2r-l) ar - 
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For large N, as r increases from m to ^ the function ^^2r){N-2r-i) mon(> 
tonically decreases from roughly jj to |. Consequently, 

i y c f 



^ 1 -{(N-m)...(N-2m)J ^ 



m<r<7V/4 



where C2 = (g) > 1. Now 



E V< E V. 

m<r<N/4 ^ m<r<oo 2 

which converges by the Cauchy criterion and, consequently, we obtain the 
desired estimate Si = 0(N~ mt ). 
To analyze £2, we first note that 

r 42 n A _(17iV/16-Ai)!_ ifA> iv 

(421) f "(JV-A0W16 + 16)! lfAl -¥' 

The proof of this estimate is similar to the proof of the estimate (4.17). 
The product of hook numbers outside the first row is at most (N — Ai)!. 
Suppose there is a total of v rows. Then 

IN 



A 2 H h A„ < 



A 2 - — + • ■ • + A15 - — < 0, 



Consequently, 



16/ V 16. 

and therefore the rows below the 16th one are of size at most That 
means that the product of hook numbers in the ^ leftmost boxes of the 
first row is no greater than + 16)!. The remaining hook lengths in the 
first row constitute Ai — distinct numbers less than or equal to N; their 
product therefore does not exceed (jv+jv/16-Ai)! 1 Application of the hook 
length formula (4.9) then completes the proof of the estimate (4.21). 

Now as Ai increases from ^ to ^f, the expression on the right-hand side 
of (4.21) decreases; therefore, in the whole range of the sum £2 we have 

A (17iV/16-3iV/4)! (5JV/16)! 

1 - (JV-3JV/4)!(JV/16 + 16)! (JV/4)!(JV/16)! ' 

and the latter expression is greater than 5 N / 16 for sufficiently large N. Con- 
sequently, 



S 2 < p(N)5~ tN / 1G < cf<Cz N = 0(c 
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for C3 = 5 4 / 16 and some C4 > 1 . 

Turning now to the sum S3, we note that 

/4\ N N 
(4.22) f x > - if Ai < Ai < — . 



Indeed, for any box a in A we have h(u) < X\ + X[ < Consequently, 
using the hook length formula (4.9) we obtain 

A N\ {N/e) N (A^ N 



f > (N/4) N > (N/4) N > U 
proving the estimate (4.22). Consequently, 

£ 3 <cl^ = 0(c 6 -") 

for C5 = (-)* and some cq > 1. 

This completes the proof of Proposition 4.2. □ 

Now Theorem 4.1 is proved by combining inequality (4.6) with Lemma 4.1, 
Theorem 4.4 and Proposition 4.2 with m = 4, together with the values of the 
remaining characters and dimensions, which are computed using Theorem 
4.2 and Theorem 4.3 and summarized in Table 1. We denote by A the set 
of partitions in the first column of Table 1. 

\\Pk *P2~U\\ < - ^ — j- j by inequality (4.6) 

<I £ ( xW(g.) )' byLemma4 , 



Table 1 
Character and dimension values 

A f x \X X (C 2 )\ \x X (Cs)\ |X A (C 4 )| |x A (C fc )|fc>5 



(TV -1,1) N-l 1 1 1 1 

(JV-2,2) H^za * 1 1 

(TV -2, 1,1) (JV - 1) 2 (JV " 2) f + 1 1 1 1 

(TV -3,2,1) N(N-2KN-4) w +1 o ! 

(TV -3, 1,1,1) (^-D(iv-2)(iv-3) w +1 1 1 

(xr_*o\ (N)(N-i)(N-5) n JV , -1 n f0, if ft = 5; 
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= i ^ ( x\c k ) x \c 2 ) \ 2 | ^/ x \ Ck ) X \C 2 ) \ 2 

Ai,A' 1 <A r -4 AeA 
= 0(iV 3/2-lA ) £ (/A) _ (1 _ 2/fc)+0(iv _ 2) 

Ai,A' 1 <A r -4 

by Theorem 4.4 and Table 1 
= 0{N^ /2+7/k ) + 0(iV^ 2 ) by Proposition 4.2; 

the first term in the last expression is at most 0(N~ 1 ^ e ), since k > 3. This 
completes the proof of Theorem 4.1. 



5. Concluding remarks. 



5.1. Extensive information is available on the distribution of cycles in 
random permutations and on Poisson-Dirichlet distribution [3, 42] and we 
will pursue the exhaustive exploitation of the consequences of Theorem 4.1 
in a subsequent publication. Here we note just two immediate corollaries 
which substantially refine results in [14, 28, 43]. 



Corollary 5.1. Let l(n) denote the number of LHT paths in a random 
cubic graph on n vertices with random orientation. Then, as n— > oo, 

(5.1) E(i(n))=log(3n) +7 + 0(1), 

(5.2) Var(/(n)) =log(3n) + 7-tt 2 /6 + o(1), 

where 7 = 0.5772 . . . is Euler's constant. Further, ; (")~|2£ n converges to stan- 

V logn 

dard normal distribution JV"(0, 1). From (2.1) we obtain that the genus is 
distributed as 1 + § — AA(logra, y/logn). 



Corollary 5.2. Let L(n) be the length of the largest LHT path. Then 

(5.3) \^e(^^\=J exp(-x-J {e- y /y)dy S j dx ~ 0.6243. 

Consequently, the expected area of the largest embedded ball in S G (T,0) 
converges to of the total surface area. 



The expression on the right-hand side of (5.3) is due to Shepp and Lloyd 
[49], who also computed the limiting distribution of L. 
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5.2. Recent numerical experiments of Novikoff [39] present convincing 
evidence in favor of the following conjecture (see [27] for a discussion of 
related conjectures and numerical results): 

Conjecture 5.1 ([39]). The distribution of the second largest eigen- 
value of the adjacency matrix of a random regular graph, suitably rescaled, 
follows Tracy- Widom GOE distribution. 

Tracy and Widom [53, 54] computed the limiting distribution function for 
the largest eigenvalue in the classical Gaussian ensembles; these distribution 
functions are expressed in terms of a certain Painleve II function and are now 
believed to describe new universal limit laws for a wide variety of processes 
arising in mathematical physics and interacting particle systems [55]. 

A dramatic consequence of Conjecture 5.1 would be that the probability of 
a random regular graph being Ramanujan approaches 0.52 as the size of the 
graph tends to infinity, corresponding to the skewness in the Tracy-Widom 
GOE distribution. 

To approach Conjecture 5.1 following the method Sinai and Soshnikov 
[50] and Soshnikov [51] in his breakthrough proof of the universality at the 
edge of the spectrum in Wigner matrices, one needs precise information for 
the number of closed walks of size up to n 2 / 3 , where n is the size of the 
graph. We hope that the results of this paper will be useful in pursuing such 
an approach. 

Acknowledgments. I started working on this problem while visiting Bob 
Brooks at the Technion in August 2001; I would like to thank the Depart- 
ment of Mathematics at Technion for its hospitality. The interplay between 
spectral geometry of surfaces and graphs, particularly the issue of the first 
eigenvalue, has been a recurrent theme in Bob's work; his papers, written 
with deep sympathy for the reader, and his continuing support and encour- 
agement have been very important to me. The results presented in this paper 
were obtained in the spring of 2002 and I communicated them to Bob a few 
months before his untimely death in September 2002. I lectured on them 
in June 2002 at the spectral geometry conference in Lexington, Kentucky 
and at the Brooks Memorial Conference in Technion in December 2003; I 
would like to thank the organizers of these conferences. After my talk at the 
Brooks Memorial Conference, Alex Lubotzky told me that Martin Liebeck 
and Aner Shalev independently obtained a result similar to Proposition 4.2; 
I would like to thank Aner Shalev for sending me their preprint [33]. I am 
grateful to Persi Diaconis, Eran Makover, Nicholas Pippenger, Anatoly Ver- 
shik, Ofer Zeitouni, Balint Virag and Martin Zerner for interest in this work 
and for stimulating discussions. 



POISSON-DIRICHLET DISTRIBUTION FOR RANDOM BELYI SURFACES 21 



REFERENCES 

[1] Alon, N. (1986). Eigenvalues and expanders. Combinatorial 6 83-96. MR0875835 
[2] Apostol, T. M. (1976). Introduction to Analytic Number Theory. Springer, New 
York. MR0434929 

[3] Arratia, R., Barbour, A. D. and Tavare, S. (2003). Logarithmic Combinatorial 

Structrues: A Probabilistic Approach. EMS, Zurich. MR2032426 
[4] Belyi, G. V. (1979). Galois extensions of a maximal cyclotomic field. Izv. Akad. 

Nauk SSSR Ser. Mat. 43 267-276. MR0534593 
[5] Bollobas, B. (1980). A probabilistic proof of the asymptotic formula for the number 

of labelled regular graphs. European J. Combin. 1 311-316. MR0595929 
[6] Bollobas, B. (1985). Random Graphs. Academic Press Inc., London. MR0809996 
[7] Bollobas, B. (1988). The isoperimetric number of random regular graphs. European 

J. Combin. 9 241-244. MR0947025 
[8] Brooks, R. (1986). The spectral geometry of a tower of coverings. J. Differential 

Geom. 23 97-107. MR0840402 
[9] Brooks, R. (1988). Some remarks on volume and diameter of Riemannian manifolds. 

J. Differential Geom. 27 81-86. MR0918458 
[10] Brooks, R. (1999). Platonic surfaces. Comment. Math. Helv. 74 156-170. 

MR1677565 

[11] Brooks, R. (2000). Some geometric aspects of the work of Lars Ahlfors. IMCP 14 
31-39. MR1786559 

[12] Brooks, R. and Makover, E. (2001). Riemann surfaces with large first eigenvalue. 

J. Anal. Math. 83 243-258. MR1828493 
[13] Brooks, R. and Makover, E. (2001). Belyi surfaces. IMCP 15 37-46. MR1890529 
[14] Brooks, R. and Makover, E. (2004). Random construction of Riemann surfaces. 

J. Differential Geom. 68 121-157. MR2152911 
[15] Buser, P. (1978). Cubic graphs and the first eigenvalue of a Riemann surface. Math. 

Z. 162 87-99. MR0505920 
[16] Buser, P. (1980). On Cheeger's inequality Ai > h 2 /A. Proc. Symp. Pure. Math. 36 

29-77. MR0573428 

[17] Buser, P. (1984). On the bipartition of graphs. Discrete Appl. Math. 9 105-109. 
MR0754431 

[18] Buser, P. (1992). Geometry and Spectra of Compact Riemann Surfaces. Birkhauser, 
Boston. MR1183224 

[19] DlACONlS, P. (1988). Group Representations in Probability and Statistics. IMS, Hay- 
ward, CA. MR0964069 
[20] DlACONlS, P. (2000). Random walks on groups: Characters and geometry. In Groups 

St. Andrews 2001 in Oxford 1 120-142. Cambridge Univ. Press. MR2051523 
[21] Diaconis, P., Mayer-Wolf, E., Zeitouni, O. and Zerner, M. (2004). The 

Poisson-Dirichlet law is the unique invariant distribution for uniform split-merge 

transformations. Ann. Probab. 32 915-938. MR2044670 
[22] DlACONlS, P. and Shahshahani, M. (1981). Generating a random permutation with 

random transpositions. Z. Wahrsch Verw. Gebiete 57 159-179. MR0626813 
[23] Fomin, S. V. and Lulov, N. (1997). On the number of rim hook tableux. J. Math. 

Sci. 87 4118-4123. MR1374321 
[24] Friedman, J. (1991). On the second eigenvalue and random walk in random d-regular 

graphs. Combinatorica 11 331-362. MR1137767 
[25] Friedman, J. (2006). A proof of Alon's second eigenvalue conjecture. Memoirs of 

the A. M.S. To appear. 



22 



A. GAMBURD 



[26] Gamburd, A. (2002). Spectral gap for infinite index "congruence" subgroups of 

SL 2 (Z). Israel J. Math. 127 157-200. MR1900698 
[27] Gamburd, A. (2004). Expander graphs, random matrices, and quantum chaos. In 

Random Walks and Geometry (V. A. Kaimanovich, ed.) 109-141. de Gruyter, 

Berlin. MR2087781 

[28] Gamburd, A. and Makover, E. (2002). On the genus of a random Riemann surface. 

Contemp. Math. 311 133-140. MR1940168 
[29] Huxley, M. (1986). Eigenvalues of congruence subgroups. Contemp. Math. 53 341- 

349. 

[30] James, G. and Kerber, A. (1981). The Representation Theory of the Symmetric 

Group. Addison-Wesley, Reading, MA. MR0644144 
[31] Jones, G. and Singerman, S. (1996). Belyi functions, hypermaps and Galois groups. 

Bull. London Math. Soc. 28 561-590. MR1405488 
[32] Kim, H. and Sarnak, P. (2003). Refined estimates towards the Ramanujan and 

Selberg conjectures. J. Amer. Statist. Assoc. 16 139-183. MR1937203 
[33] Liebeck, M. E. and Shalev, A. (2004). Fuchsian groups, coverings of Riemann 

surfaces, subgroups growth, random quotients and random walks. J. Algebra 

276 552-601. MR2058457 
[34] Lubotzky, A. (1994). Discrete Groups, Expanding Graphs and Invariant Measures. 

Birkhauser, Basel. MR1308046 
[35] Lubotzky, A., Phillips, R. and Sarnak, P. (1988). Ramanujan graphs. Combi- 

natorica 8 261-277. MR0963118 
[36] Luo, W., Rudnick, Z. and Sarnak, P. (1995). On Selberg's eigenvalue conjecture. 

Geom. Fund. Anal. 5 387-401. MR1334872 
[37] McKay, B. (1981). The expected eigenvalue distribution of a large regular graph. 

Linear Algebra Appl. 40 203-216. MR0629617 
[38] Mulase, M. and Penkava, M. (1998). Ribbon graphs, quadratic differentials on 

Riemann surfaces, and algebraic curves defined over Q. Asian J. Math. 2 875- 

919. MR1734132 

[39] Novikoff, T. Asymptotic behavior of the random 3-regular bipartite graph. 
Preprint. 

[40] Okounkov, A. (2002). Symmetric functions and random partitions. In Symmet- 
ric Functions 2001: Surveys of Developments and Perspectives (S. Fomin, ed.). 
Kluwer, Dordrecht. MR2059364 

[41] Pinsker, M. (1973). On the complexity of concentrator. 7th Annual Teletrafic Con- 
ference, 318/1-318/4, Stockholm. 

[42] Pitman, J. Combinatorial stochastic processes. Technical Re- 
port 621, Dept. Statistics, U.C. Berkeley. Available at 
http : //stat-www.berkeley . edu/users/pitman/bibliog.html . 

[43] Pippenger, N. and Schleich, K. (2006). Topological characteristics of random tri- 
angulated surfaces. Random Structures Algorithms 28 247-288. MR2213112 

[44] Randol, B. (1974). Small eigenvalues of the Laplace operator on compact Riemann 
surfaces. Bull. Amer. Math. Soc. 80 996-1008. MR0400316 

[45] Regev, A. (1981). Asymptotic values for degrees associated with strips of Young 
diagrams. Adv. in Math. 41 115-136. MR0625890 

[46] Sarnak, P. (2004). What is ... an expander? Notices Amer. Math. Soc. 51 762-763. 
MR2072849 

[47] Sarnak, P. and Xue, X. (1991). Bounds for multiplicities of automorphic represen- 
tations. Duke Math. J. 64 207-227. MR1131400 



POISSON-DIRICHLET DISTRIBUTION FOR RANDOM BELYI SURFACES 23 



[48] Selberg, A. (1965). On the estimation of Fourier coefficients of modular forms. Proc. 
Symp. Pure Math. 8 1-15. MR0182610 

[49] Shepp, L. A. and Lloyd, S. P. (1966). Ordered cycle lengths in a random permu- 
tation. Trans. Amer. Math. Soc. 121 340-357. MR0195117 

[50] Sinai, Ya. and Soshnikov, A. (1998). A refinement of Wigner's semicircle law in 
a neighborhood of the spectrum edge for random symmetric matrices. Fund. 
Anal. Appl. 32 114-131. MR1647832 

[51] Soshnikov, A. (1999). Universality at the edge of the spectrum in Wigner random 
matrices. Comm. Math. Phys. 207 697-733. MR1727234 

[52] Tsilevich, N. V. (2000). Stationary random partitions of a natural series. Theory 
Probab. Appl. 44 60-74. MR1751188 

[53] Tracy, C. A. and Widom, H. (1994). Level-spacing distribution and the Airy kernel. 
Comm. Math. Phys. 159 151-174. MR1257246 

[54] Tracy, C. A. and Widom, H. (1996). On orthogonal and symplectic matrix ensem- 
bles. Comm. Math. Phys. 177 727-754. MR1385083 

[55] Tracy, C. A. and Widom, H. (2002). Distribution functions for largest eigenvalues 
and their applications. In Proceedings of the International Congress of Mathe- 
maticians, Beijing I (L. I. Tatsien, ed.) 587-596. Higher Education Press, Bei- 
jing. MR1989209 

[56] Vershik, A. M. (1989). A partition function connected with Young diagrams. J. 
Soviet Math. 47 2379-2386. MR0947326 

[57] Vershik, A. M. and Shmidt, A. A. (1977). Limit measures that arise in the asymp- 
totic theory of symmetric groups. I. Teor. Verojatnost. i Primenen. 22 72-88. 
MR0448476 

[58] Vershik, A. M. and Shmidt, A. A. (1978). Limit measures that arise in the asymp- 
totic theory of symmetric groups. II. Teor. Verojatnost. i Primenen. 23 42-54. 
MR0483019 

[59] Wormald, N. C. (1999). Models of random regular graphs. In Survey in Combina- 
torics 1999 239-298. Cambridge Univ. Press. MR1725006 

Department of Mathematics 

University of California 

1156 High Street 

Santa Cruz, California 95064 

USA 

E-MAIL: agamburdOucsc.edu 



